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We study the dynamics of transverse oscillations of a suspended carbon nanotube into which 
current is injected from the tip of a scanning tunneling microscope (STM). In this case the corre- 
lations between the displacement of the nanotube and its charge state, determined by the position- 
dependent electron tunneling rate, can lead to a "shuttle-like" instability for the transverse vibra- 
tional modes. We find that selective excitation of a specific mode can be achieved by an accurate 
positioning of the STM tip. This result suggests a feasible way to control the dynamics of this 
nano-electromechanical system (NEMS) based on the "shuttle instability" . 

PACS numbers; 85.35.Kt, 85.85.-|-j 



There are several reasons for the considerable cur- 
rent interest in nano-electromechanical systems (NEMS), 
both for technological applications and fundamental re- 
search. The peculiar combination of several features 
such as high vibrational frequencies and small masses 
which characterize most NEMS makes these systems 
very suitable for the realization of new measurement 
tools with extremely high sensitivity in mass sensing and 
force microscopy applications [III4I. Furthermore, the 
mechanical elements of the NEMS (typically cantilevers 
or beams) are considered the most promising structures 
where quantum features of motion could be experimen- 
tally detected [1]. 

The physical basis for many of the interesting function- 
alities of NEMS is the strong interplay between mechani- 
cal and electronic degrees of freedom In the par- 
ticular case of a nano-electromechanical single-electron 
transistor device having a metallic dot as movable part, 
the equilibrium position of the dot can become unstable 
as a consequence of the electromechanical coupling. In 
this case the dominant mechanism for the transport of 
charge is based on the oscillations of the dot which can 
"shuttle" the tunneling electrons across the system [^Q. 

The typical set-up for many NEMS includes a spatially 
extended movable element such as a suspended carbon 
nanotube, whose dynamics has been demonstrated to be 
characterized by a number of different vibrational modes 
[Toj . The relevance of many mechanical modes in the 
transport of charge suggests that the variety of effects 
due to the electromechanical coupling in suspended car- 
bon nanotube-based NEMS may be even richer than in 
the ordinary "shuttle" system. Jonsson et al. have shown 
TllfT^flst that if extra charge is injected into the movable 
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part of the device from the tip of a scanning tunneling mi- 
croscope (STM) a nano-electromechanical "shuttle-like" 
instability can be induced for the transverse vibrational 
modes of the nanotube. 

The selective promotion of the electromechanical insta- 
bility for different vibrational modes provides an interest- 
ing perspective for probing the dynamics of NEMS. Here 
we show that such selective excitation can be achieved 
by means of local injection of electric charge. The main 
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Figure 1: Sketch of the model system considered. The dis- 
tance of the suspended carbon nanotube from the STM tip 
affects the electron tunneling rate between them, while the 
tunneling rate between the nanotube and the leads is constant 
(the vibration amplitude and the STM-nanotube equilibrium 
distance d shown are exaggerated for clarity). 



idea is to optimize the electromechanical coupling for the 
mode(s) which we want to make unstable. The local char- 
acter of the electric charge injection makes the selective 
excitation of the nanotube transverse modes possible by 
varying the position of the STM tip. 

We will consider here the same device analyzed by 
Jonsson et al. since it provides a convenient set-up to 
control the electromechanical coupling of different vibra- 
tional modes. The system is sketched in Fig. ^ and 
it consists of a suspended metallic carbon nanotube in 
tunneling contact with an STM tip and one supporting 
lead. 

We take the z-axis along the nanotube axis, while its 
cross section lies in the xy plane. The STM tip is put 
over point (0, 0, zq) and its distance from the nanotube 
at equilibrium is c? ~1 nm. 

In order to describe the motion of the nanotube we 
model it as a classical elastic beam of length L clamped 
at both ends and focus on its flexural vibrations. 

The motion of the nanotube in the xz plane can be 
described through u{z,t), its displacement along the x 
axis from the static equilibrium configuration. If the am- 
plitude of the oscillations is small enough for linear elas- 
ticity theory to be valid, the time evolution of u{z,t) is 
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determined by the following equation \l£ 
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In Eq. IT]), p is the carbon nanotube density, S the cross 
section, E the Young modulus, / the cross section mo- 
ment of inertia, N{t) the number of extra electrons on the 
nanotube at time t and F^^ the x component of the exter- 
nal force (per unit length) generated by the electrostatic 
interaction between the STM tip and the nanotube. 

The precise spatial distribution of depends on the 
details of the geometric structure of the tip apex. How- 
ever, a simple electrostatic analysis indicates that for 
\z — zo\ ^ R, d (where R < 10~^ m is the effective 
size of the STM tip), the force decays at least as 
\z — zo|~'^- Therefore the influence of the metallic leads 
can be ignored as long as the STM tip is not too close to 
them and we can write F°'(iV, z, zq) ~ F°^{N, z — zq). 

The displacement field u{z,t) and the force per unit 
length can be expressed as linear combinations of 
eigenf unctions ipj{z) of the operator d^/dz^ with the 
boundary conditions u(0,f) = u{L,t) — 0, du/dz{0,t) = 
du/dz{L,t) = 0. 

The expansion of u{z, t) and over the complete set 
of functions ^Pj{z) makes it possible to decompose Eq. 
(fT|) into a set of equations of motion for the eigenmode 
amplitudes Xj{t), which can be written in hamiltonian 
form by introducing the conjugate momenta T^j{t)'. 
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+ Ij^j + rriLUjXj = Lff{N{t), {xj}, zq). 



(2a) 
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In Eqs. fj^ are the coefficients in the expan- 

sion of F^^ over the complete set of functions (pj{z) 
(which are chosen to be dimensionless), m — pSL is 
the mass of the nanotube and ujj are the frequencies 
of the transverse vibrational modes, given by luj = 
^ykjEI/ {pSL"^), where the eigenvalues kj can be found 
by solving: cosfcj^^coshfcj^^ — 1. 

We introduced in Eq. (j2bp a phenomenological damp- 
ing force for each mode, — 7j7rj, where 7j has the di- 
mension of inverse time. The motion of the nanotube is 
inevitably affected by dissipative mechanisms, which can 
be related to its coupling to the environment and to sev- 
eral internal processes. We will later consider a general 
form for jj which can be used to describe the damping 
induced by a wide class of phenomena. 

An approximate expression for the force coefficients ff 



in Eq. pb|) can be found through some physical consider- 
ations on the characteristic lengths of the system. Since 
the eigenfunctions fj{z) vary appreciably only over dis- 
tances of the order of L, we can express each /j' as a sum 
of a sharply localized contribution at zq plus a correction: 



Lff = F^\N{t), z' - zo)ip,{z')dz' = 

J-L/2 

= eNit)£ipjizo) + 0{R^/L^), (3) 

where f is a phenomenological parameter which provides 
the magnitude of the effective electrostatic field between 
the STM and the nanotube. 

The size of the correction in Eq. ([3]) can be estimated 
in terms of the characteristic lengths of the system: for 
the typical values L ~ 10~^ m and R ~ 10~^ m, the con- 
dition R^ / ^ 1 is fulfilled and that defines the range 
of validity of the approximation Lf^^ ~ eN{t)£(pj{zQ) 
which we will use from now on. 

For what concerns the transport of charge the system is 
equivalent to a double tunnel junction, having one poten- 
tial barrier localized between the STM tip and the nan- 
otube and the another one between the nanotube and 
the leads. In our analysis we will consider the case of 
electrons for which the inverse characteristic time of de- 
coherence is much shorter than the tunneling rates, so 
that the description of tunneling as a stochastic (rather 
than coherent) process is sufficient. 

Following the approach presented in for a point- 
like oscillator coupled to a single-electron transistor, we 
define a probability density in the phase space of the 
system Pn{x, tt, t) such that Pn{x, tt, t)dxid'Kidx2d'K2 ■ ■ ■ 
(with all the Xj and nj scaled by proper dimensional fac- 
tors) is the joint probability that at time t there are N 
electrons in excess on the nanotube while the eigenmode 
amplitudes and momenta {xj}, {tTj} take values in the 
phase space region defined by Ylj dxjd-Kj = dxdir. 

We consider the system in the Coulomb blockade 
regime and limit to one the maximum number of extra 
electrons on the nanotube, therefore only the probabil- 
ity densities Po{x,7T,t) and Pi{x,7T,t) play a role in the 
description of the nanotube dynamics. 

The coupling between the mechanical and electronic 
degrees of freedom arises because the tunneling rate be- 
tween the STM tip and the nanotube, rsTM-NT(a?7 ^o) 
is affected by their relative distance at point zq: 
rsTM-CNT(^, 2;o) = roexp((-c; + Y.j^jV]{za))/>^), 
where A is the effective tunneling length of the STM- 
nanotube junction and Fq is a constant. The factor 
rg(ci) = roexp(— d/A) is the tunneling rate that would 
characterize the junction if the motion of the nanotube 
could be neglected. The tunneling rate between the nan- 
otube and the leads, Fnt-l does not depend on the nan- 
otube displacement. 

We remark that all the tunneling rates are generally 
functions of the bias voltage. However, since here we al- 
ways assume V fixed at some value we never explicitly 
indicate the dependence on V . In order to be consis- 
tent with the condition of single-electron charging of the 
nanotube, the electron temperature and the bias voltage 
have to be in the range: k^T < Eq < eV < 2Ec, where 
Eq is the energy required to add one extra electron to 
the nanotube. 
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The time evolution of the probabiHty densities 

P+ (x, TT , t) = Pi (x, TT , t) + Pq (x, TT , t) and P_ (x, TT,t) = 

Pi {x, TT, t) — Po{x, TT , t) is determined by the equations: 



dP+ 

dP- 
~dt 



{Ci+C2){P+ 
{Ci+C2){P+ 



P-) = 



P-) = 



r_(x; zo)P+ + r+(f; zq)P^, 



(4a) 



(4b) 



where T+{x; zq) = rsTM-NT(x; zo) + rNT-L, r_(f; 2:0) = 
rsTM-NT(x; ^o) — Tnt-l and the Liouvilhan operators 
Ci and £2 are defined as foUows: 
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m dxi 
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From Eqs. ([4]) we can derive the equations of motion 
for any dynamical variable averaged over the probability 
densities P+ and P_: ((...))„ = f{...)Pa{x,Tf,t)dxdTr, 
where a — ±. The set of equations of motion for the first 
moments =p-(t), {xj)a, {'^j)a is: 



dt m 

d{nj)a 



dt 

dp^ 

~dr 



(5a) 

Hi 

-muj'^{xj)a - l]{'^i)a + e£ipj{zo)pi (5b) 



(r_(f;zo))+ - (r+(f; zo))-, 



(5c) 



The set of equations (O is not closed because the ex- 
ponential form of the tunneling rate rgTM-NT(x; •^o) in- 
troduces a coupling between the first and all the other 
moments. However in the limit of small oscillation am- 
plitudes we can expand rsTM-NT(2;; -Zo) to first order in 
Xj/X and that reduces ([5]) to a closed set of linear equa- 
tions. 

The static solution of the linearized equations of mo- 
tion, {Trj)a = 0, {x°')a = x'j, p- = p, where and 
p are constant, describes the nanotube as a slightly bent 
beam at rest. The stability of this solution can be investi- 
gated by substituting the expressions {£,k)+ = Cfe+^fc^'^''* 
(where is any of the dynamical variables (x^), (TTfe), p_ 
and Afc is constant) in the linearized equations of motion 
and solving for the exponents Pk 

This procedure leads to an algebraic equation which 
in general cannot be solved analytically. However, if 
the dimensionless parameters £k = e£(p1{zo) / {muj'^X) are 
much smaller than 1 and the nanotube is only weakly 
damped (7^ ^ LUk), we can look for exponents of the 
form iwfc -1- Sk, with \dk\ (jJk and derive an ana- 

lytical expression for the Sk which up to the first order 
in all £k and jk reads: 
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rs(c^)rNT-L w^efe(zo) 



2r;(rf) 



Tfid) 



l + i 



r;(rf) 

2ujk 



(6) 

where r\,{d) = T^{d) + Tnt-l- The condition e/c < 1 
can be taken as a definition of the weak electromechanical 
regime, since it implies that the shift in the equilibrium 
position of the nanotube at point zq when it is charged 
by one extra electron, Su = e£ J2 j (^0)/ (tow^), is much 
smaller than the tunneling length A. For realistic values 
of the parameters which are consistent with the condi- 
tions of validity of our model {V ^ 0.1 V, the capac- 
itances of the tunnel junctions Cstm-nt ^ C'nt-l ~ 



10- 



-18 



F as reported in [17], A 



10- 



-10 



m, Uk 



MHz), m ^ 10^^^ kg, L ^ 10^^ m, the regime of weak 
electromechanical coupling is attained: £k{zo) ^0.1 for 
all the modes at every position zq- 

The sign of the real part of 5k in Eq. ^ determines 
the stability of the static solution for the fc— th average 
mode amplitude. If 5He[(5fc] > then {xk)+ increases ex- 
ponentially in time, hence the static solution for the fc-th 
mode is unstable. This is the signature of a "shuttle- 
like" electromechanical instability. On the other hand, if 
9^e[(5fc] < the energy pumped into the vibrational mode 
by the electrostatic field is not able to compensate the 
loss due to dissipation and after a time interval of the or- 
der of 1 /7fe the /c-th mode amplitude decays to its static 
value. 

For fixed values of 7^, £ and the sign of 9^e[(5fe] 
becomes a function of zq and rg(d), i.e. it depends only 
on the position of the STM tip in the xz plane. The set 
of values of and Fq [d) for which the real part of 5k is 
positive defines the instability region for the fc-th mode 
amplitude in the plane (zo,rg(d)). 

In order to map the instability regions we have to spec- 
ify an analytic expression for the damping rates 7fc . The 
dissipation of mechanical energy in NEMS can take place 
through several mechanisms. However, in spite of the 
variety of the dissipative processes in solids, their effect 
on the NEMS performance can be described by consid- 
ering the retardation induced in the NEMS response to 
mechanical perturbations (which adds to the "instanta- 
neous" elastic behaviour). 

In order to include this effect in our model we follow 
the approach introduced by Zener and formally replace 
the Young modulus with a frequency dependent complex 
function which results in the following expression for the 
damping rates: 7/c = Yuj1,t/{1 + ojf.T'^) |14 |. 

A large class of dissipative phenomena in solids (e.g. 
thermoelasticity, dislocations and defects dynamics) can 
be parametrized though the dimensionless coefficient Y 
and the relaxation time r, which both depend on tem- 
perature as well as on several geometric and material 
properties p^ . 

We first consider the limit in which the characteristic 
inverse time of the retarded mechanical response is much 
smaller than the frequencies of the nanotube eigenmodes: 
\/t <^ ujk- In this case the damping term in Eq. ^ does 
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Figure 2: (Color online) Regions of instability in the param- 
eter plane {zq, ro(ci)) for the first three bending modes when 
the inverse characteristic time of the nanotube mechanical re- 
sponse to the external force is smaller than the frequencies of 
the nanotube vibrational modes, 1/r = O.ltJi. The regions 
computed from the real part of the analytic expression (|6]) 
(solid lines) are compared with the results obtained by the nu- 
merical solution of the linearized equations of motion (mark- 
ers). The parts of the plot filled with colors correspond to the 
values of (zo, ro(d)) for which a single vibrational mode is ex- 
cited. The other relevant parameters are e£ / {mXuf) — 0.1, 
Fnt-l = Swi, Y = W-\ 
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Figure 3: (Color online) Regions of instability in the param- 
eter plane (zo, ro(d)) for the first three bending modes when 
the inverse characteristic time of the nanotube mechanical 
response to the external force is larger than the frequencies 
of the nanotube vibrational modes, 1/r = lOoji. The re- 
gions computed from the real part of the analytic expression 
((B]) (solid lines) are compared with the results obtained by 
the numerical solution of the linearized equations of motion 
(markers) . The region of instability for the first mode includes 
all the values of (zo, rS(d)) for which the other modes are un- 
stable, therefore no selective excitation can be attained in this 
regime. The other relevant parameters are 
Fnt-l = 5lji,F = 10-'' . 



not depend on the frequency and the dissipation rate is 
the same for all the modes. 

In Fig. the instability regions determined by Eq. 
^ for the first three modes are plotted together with the 
real parts of the exponents /3fc obtained from the numer- 
ical analysis of the linear stability problem. The areas 
filled with colours correspond to the values of (^q, rQ((i)) 
for which only a single vibrational mode is unstable, while 
the regions where two or three modes are unstable are left 
blank. 

The physical picture presented in Fig. ^ changes 
drastically in the opposite limit, 1/t ^ uJk, as can be 
seen in Fig. ([3]). In this case the first mode is character- 
ized by the smallest dissipation rate, 71 ^ 7/c k ^ 1, 
therefore if any of the modes is unstable, also the first one 
is unstable. That excludes the possibility of promoting a 
selective instability in the limit l/r ^ tj^ . 

The dynamical behaviour of the nanotube in the 
regime of single-mode instability is qualitatively the same 
of the ordinary "shuttle" system [g]. The amplitude of 
the oscillations increases exponentially until it reaches a 
certain steady value which depends on the parameters of 
the system. This transition is characterized by a large 
enhancement of the current (with respect to the static 
tunneling regime) that can be experimentally detected 
by measuring the current flowing through the device for 
different positions of the STM tip. 

In the phase space of the system the dynamical state in 



this situation is described by a limit cycle, that is an iso- 
lated closed trajectory characterized by finite amplitude 
oscillations [l^. Here we have shown that the frequency 
of this stable self-oscillating state can be selected among 
the whole set of nanotube resonant frequencies through 
an accurate positioning of the STM tip. 

In conclusion in the present work we studied the dy- 
namics of the flexural vibrations of a suspended car- 
bon nanotube in which extra electrons are injected at 
a position-dependent rate. We showed that a localized 
constant electrostatic field can excite many transverse 
vibrational modes of the nanotube into a "shuttle-like" 
regime of charge transport. For a fixed bias voltage and 
in presence of dissipative processes with inverse charac- 
teristic times much smaller than the frequencies of the 
nanotube vibrational modes, we found that it is possi- 
ble to induce a selective instability through an accurate 
positioning of a STM. It thus seems possible to extend 
the approach followed here to other systems character- 
ized by a non trivial coupling between charge transport 
and mechanical degrees of freedom. 

The author wants to thank L. Y. Gorelik, R. I. 
Shekhter and M. Jonson for fruitful discussions and sup- 
port. Partial financial support from the Swedish VR and 
from the Faculty of Science at the University of Gothen- 
burg through its "Nanoparticle" Research Platform is 
gratefully acknowledged. 



5 



^ M. L. Roukcs, Tech. Rep., Transducer Research Founda- 
tion (2000), arXiv:cond-mat/0008187v2. 

^ M. P. Blencowe, Contemporary Physics 46, 249 (2005). 

^ K. C. Schwab and M. L. Roukes, Physics Today 58, 36 
(2005). 

" P. Poncharal, Z. L. Wang, D. Ugartc, and W. A. doHccr, 

Science 283, 1513 (1999). 
^ S. T. Purcell, P. Vincent, C. Journet, and V. T. Binh, 

Phys. Rev. Lett. 89, 276103 (2002). 
® S. Sapmaz, Y. M. Blanter, L. Gurevich, and H. S. J. van der 

Zant, Phys. Rev. B 67, 235414 (2003). 

V. Sazonova, Y. Yaish, H. Ustiincl, D. Roundy, T. A. Arias, 

and P. L. McEuen, Nature 431, 284 (2004). 
® L. Y. Gorelik, A. Isacsson, M. V. Voinova, B. Kascrno, 

R. I. Shekhter, and M. Jonson, Phys. Rev. Lett. 80, 4526 

(1998). 

^ A. V. Moskalcnko, S. N. Gordccv, O. F. Kocntjoro, P. R. 

Raithby, R. W. French, F. Marken, and S. E. Savel'ev, 

Phys. Rev. B 79, 241403(R) (2009). 
° A. K. Hiittel, M. Pott, B. Witkamp, and H. S. J. van der 



Zant, New J. Phys. 10, 095003 (2008). 
" L. M. Jonsson, L. Y. Gorelik, R. L Shekhter, and M. Jon- 
son, Nano Lett. 5, 1165 (2005). 

L. M. Jonsson, L. Y. Gorelik, R. I. Shekhter, and M. Jon- 
son, New J. Phys. 9, 90 (2007). 
^'^ L. M. Jonsson, F. Santandrca, L. Y. Gorelik, R. I. 
Shekhter, and M. Jonson, Phys. Rev. Lett. 100, 186802 
(2008). 

A. N. Cleland, Foundations of nanomechanics (Springer- 
Verlag, 2003). 

A. D. Armour, M. P. Blencowe, and Y. Zhang, Phys. Rev. 
B 69, 125313 (2004). 

S. Strogatz, Nonlinear dynamics and chaos: with ap- 
plications to physics, biology, chemistry and engineering 
(Perseus Books Group, 2001). 

B. J. LeRoy, S. G. Lemay, J. Kong, and C. Dekker, Appl. 
Phys. Lett. 84, 4280 (2004). 

A. S. Nowick and B. S. Berry, Anelastic relaxation in crys- 
talline solids (Academic Press, 1972). 



